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Abstract 

This paper describes how coherent backscattering is altered by an external magnetic field. In the theory presented, magneto- 
optical effects occur inside Mie scatterers embedded in a non-magnetic medium. Unlike previous theories based on point-like 
scatterers, the decrease of coherent backscattering is obtained in leading order of the magnetic field using rigorous Mie theory. 
This decrease is strongly enhanced in the proximity of resonances, which cause the path length of the wave inside a scatterer 
to be increased. Also presented is a novel analysis of the shape of the backscattering cone in a magnetic field. 




Introduction 

^^he enhancement of backscattering in propagation of waves in a random medium is a well documented topic. Weak localization 
CJheory explains how interference effects between direct and reverse scattering events produce Coherent Backscattering (CBS). 

The main features of CBS are insensitive to many aspects of the statistics of the inhomogeneities. Even absorption does not alter 
T-^h]e relative strength of the interference ^ . The interference is only affected by the properties of the medium with respect to the 
^Reciprocity principle as first noted by Golubentsev Q. MacKintosh and John Q analyzed how CBS is altered by Faraday rotation 
^and natural optical activity in a medium of inhomogeneities smaller than the wavelength, and in the diffusion approximation. 
Rising a method based on point-like scatterers, Van Tiggelen et al. extended these ideas and discussed the anisotropy induced by 
^ magnetic field in light diffusion . The case of non-magnetic Mie scatterers embedded in a Faraday-active medium has been 
at ddied by means of Monte Carlo simulations by Martinez et al. |^] . Maradudin et al. considered specifically the two dimensional 
Obojierent backscattering of light from a randomly rough surface in the presence of a magnetic field Experimentally, CBS in 
magnetic field has been studied by Erbacher, Lenke and Maret [Q, and some of their results will be discussed here. 
In section ^, the main results of a recent calculation of the T-matrix for a Mie scatterer in a magnetic field are presented 
~|, and serve in section ^ as the building block to study diffusion of light in a magnetic field. After having detailed the main 
eatures of the Faraday effect for multiple Rayleigh and Mie scattering, the modification of the line shape of CBS in a medium 
'^Ji^ith finite-size scatterers in a magnetic field is investigated in the last section. 

a : 
o . 

T-matrix in a magnetic field 

> . 

•^n this paper, cq — \ has been set. In a magnetic field, the refractive index is a tensor of rank two. For the standard Mie 
r^roblem, it depends on the distance to the center of the sphere r, which has a radius a, via the Heaviside function 8(r — a), that 
j^quals 1 inside the sphere and outside, 

e(B,r)-I= [(eo-l)I + £f *]e(r-a). (1) 

In this expression, Eq = m? is the value of the normal isotropic dielectric constant of the sphere of relative index of refraction 
m, and 

ep = 2mVoB/uj (2) 

is a dimensionless coupling parameter associated with the amplitude of the Faraday effect (Vq being the Verdet constant, B the 
amplitude of the magnetic field and uj the frequency). The antisymmetric hermitian tensor = ietjkBk has been introduced 
(the hat above vectors means that the vectors have been normalized). The Mie solution depends on the dimensionless size 
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parameters x — toa and y — mx. In this paper, only non-absorbing media are considered so that m and ep are real- valued. Since 
ep ~ 10~* in most experiments, a perturbational approach is valid. The part of T that is independent of the magnetic field is 
denoted T", the part of the T-matrix linear in B is T^, and the second order correction T^. 

Two important symmetry relations must be obeyed by a T-matrix of a scatterer in a magnetic field. The first one is parity 
symmetry and the second one is reciprocity j^: 

r_kcr,-k'CT' (B) — TkCT,k'<T'(B), (3) 



T-k'-.',-k-.(-B) = Tk.,kV'(B). (4) 

It is important to note that cr(— k) = — cr(k), i.e. a indicates the helicity. In particular, relations (|^, ^ must hold for T^. 
satisfies Eq. without a minus sign for B and obeys the standard reciprocity principle. For this reason, will not contribute 
to the suppression of the backscattering cone considered in this article, the next order being T'^. 

Because is linear in B, it can be constructed by considering only three special cases for the direction of B. If k and k' are 
not collinear, the unit vector B can be decomposed in the non-orthogonal but complete basis of k, k' and g = k x k'/|k x k'|. 
This results in §, 

, (B-k)(k-k')-B-k' 1 ,~ (B-k')(k-k')-B-k_, 
Tkk' = — T^^k' (B = k' + -^^ T^k' B = k) 

(k • k')2 - 1 (k • k')2 - 1 

+ (B.g)Tik'(B = g), (5) 
With respect to an helicity base, takes the form, 

r,V(B = k) = -[R^{e)a + R^ieyi (6) 

Tl,iB = k') = -[RAOy + R2{9)ct], (7) 

UJ 

with 

RlW = -^y] r7r ', ^^ [Cj^J,l (^) + ^^J^l (8) 

[Dj7rj,i{e) + CjTj,,i0)] (9) 
where the coefficients Cj and Vj have been defined in Ref . M . For T^^., (B = g) , the following expression has been found 
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T^,,{B = i) = -(Q,+aa'Q2) (10) 
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For k = k', 



Qr = -4i?r, = (11) 

27r 

TU = 1'-^i(0), (12) 



ith 

ii;i(0) = i?2(0) = -— V(2J + 1) (Cj + Vj). (13) 
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3 Diffusion of light in a magnetic field 



The equation of radiative transfer aims to describe the propagation of the average intensity in multiple scattering, but violates 
reciprocity The reason is that it does not contain the most-crossed diagrams (C^f^i) responsible for CBS. By definition, 
the four-rank tensor Lijki linearly connects the field correlations of incident and outgoing fields. This can be represented by a 
diagram that starts on the first scattering event, and which has two lines, corresponding to the propagation of the field and to 
the complex conjugate of the field as can be seen in Fig. (|l|). For L^^,; both fields propagate in the same magnetic field B as 
opposed to for which both fields propagate with opposite sign for the magnetic field B [Q. The relations deduced from the 
application of the reciprocity relation (|^) are 

C-,(p,p',q,B) . I^,, (E^, p:^,P + p',b) , (14) 

and 

L±,,(p, p', q, B) = L%, (-p', -p, -q, -B) . (15) 

On long length scales (q — > 0) and for stationary situations (il = 0), the ladder diagrams i"^,; take the following hydrodynamic 
form [0, 

r- f I TjN 27rdifc(p,q,B)rfij(-p',-q,-B) 

i,^.,,(p,p ,q,B) = ^ ■ (16) 

The symmetry of the numerator is imposed by the reciprocity principle; A is a scalar dimensionless parameter that elucidates 
the breaking of reciprocity by the magnetic field in multiple scattering. In Eq. (|l^) , D is the conventional diffusion coefficient 
for radiative transfer, and (. is the scattering mean free path, which is a typical distance between two scattering events. This 
diffusion constant was shown theoretically and experimentally to depend on the square of the magnetic field. These 
corrections are not considered here since the dephasing parameter X/(D£) is only discussed to second order in the magnetic field. 
The enhancement factor of coherent backscattering is a ratio of a coherent contribution, described by the most-crossed diagrams, 
over an incoherent contribution represented by the ladder diagrams. The parameter A plays the role of an "absorption" term for 
the coherent contribution only, and is therefore responsible for the decrease of the coherent backscattering cone in a magnetic 
field. This parameter A can be expressed as the square of a Faraday dephasing angle, which is the product of a Verdet constant, 
the magnetic field strength, and a characteristic length scale for the propagation of light in the medium. Since the scattering 
mean free path £ is a natural and experimentally relevant length scale for coherent backscattering, it is possible to define an 
effective Verdet constant Ve/f from a relation derived in , 



A = ^V^f^B^e'^. (17) 



Equation ([17|) was obtained by MacKintosh and John |^ , who considered a situation where the scatterers and the outside medium 
both have a Verdet constant Vef / . 

To first order in the magnetic field, the tensor dik defined in Eq. (Um must take the form 



d(p,q = 0,B) =I + d2$. (18) 

where d2 describes persisting polarization effects in diffuse scattering due to the magnetic field ||T^. For low density n of particles, 
^2 can be determined from the Bethe-Salpether equation 

Lpp,(B) - nTpp,(B)T;p,(-B) +n^Tpp„(B) • G(p",B) • G*(p", -B) • T^p^-B) • L-„^,(B), (19) 

p" 

where G(p, B) denotes the ensemble averaged Dyson Green's tensor (see Appendix A for details of notations), and the asterisk 
denotes hermitian conjugation in polarization space. This equation is slightly more complicated than the one for L^p' |pT| 
because it involves both B and — B due to the definition (p^. Inserting Eqs. (|l^) and ( [l^ ) into Eq. ([l^) and expanding to first 
and second order in the magnetic field using Eq. (|^), fixes d2 and A rigorously. The final result reads 

^2^»( ^ ~v], (20) 

\ 1 - (cost'),, / 
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(21) 



(22) 



(23) 



•^'^Q scatt J —I 

Equations ( pO[-p^ ) contain the main results of this paper and wiU be discussed in sections |^ to |[ The factor 1/(1 — (cos 9)^) in 
Eqs. ( |20| ) and ( pl| ) is related to a depolarization length which will be introduced in the next section. In equation (pT]), A contains 
two terms: the first one is seen to dominate for finite size particles {x ^ 1) whereas the second one prevails for small particles 
(x 1). Sections |^ and ^ will discuss these two contributions respectively. Far from resonances, the first term corresponds to 
the Faraday rotation of the wave inside the scatterer, whereas the second represents the Faraday rotation between two scattering 
events. Although the medium outside the scatterers is not Faraday active, the Faraday rotation from one scatterer to the next 
can be defined in the framework of an effective medium theory using the parameter ry, which will be defined in section 



4 Depolarization length 

In this section, the factor 1/(1 — (cos^)^), which tends to amplify the second term of Eq. is defined. The term A2 is easily 
seen to dominate for small particles. The total cross-section of one particle is given by ll^. 



hcaU = \ [ i[Si\' + \S2\')dfi, 



(24) 



with /i — cos 6*. The transport mean free path is defined as £* = £/{l — (cosd)), where the asymmetry parameter {cos6) is given 

by 



1 



{[s,\^ + \S2\'')^ld^l. 



C scatt J — I 

Likewise, a depolarization length can be defined as £dep = ^/(l ^ (cos^)^), where 
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"P x^Q 



scatt J ^1 

Rayleigh scattering has a forward-backward symmetry so that {cos6) — 0. However, scattering is not isotropic due to the 
polarization and one can easily show that for Rayleigh scatterers {cosO)^ — 0.5. In the limit of large forward scattering and for 
m = 1.33, both asymmetry parameters, (cos 6*) and {cosO)^ tend towards a limit close to 0.85, as shown in Fig. ^t). In the 
forward direction the differences between the two states of polarization vanish since Si (0) = S2 (0) . 

The well-known oscillations and ripple structure jl3| of the asymmetry parameter (cos 6*) are also present in {cos 6) p. As 
shown in Fig. (||b), for a relative high value of the relative index of refraction (m = 2.73 corresponding to Ti02 ), the 
asymmetry parameters {cos 9) and (cos^)^ may take negative values, which can be seen in this particular case near x — 2. In 
this very particular case, where the scattering is essentially in the backward direction, the characteristic length for the loss of 
the polarization Ij^gp can be smaller than the characteristic length for the loss of the phase in multiple light scattering, which is 



5 Faraday rotation for multiple Rayleigh scattering 

Equation (|2^) for d2 and the second term in Eq. ( |2l| ) can be understood from an "effective medium" theory, valid for Rayleigh 
scatterers. The real part of the forward scattering amplitude (|l2|) is associated with the Faraday effect and the imaginary part 
with magneto-dichroism (i.e different absorption for different circular polarization) of an ensemble of Faraday-active scatterers. 
For a dilute system, the antisymmetric part of this effective refractive index Sa is defined as 
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27m , , , ^ 

ea = . 25 

For the real part of this effective refractive index, the dimensionless parameter 77 is defined by 

f] = 3?e(£a)H. (26) 

Equations ( |20| ) and ( p3| ) involve this parameter rj, which represents a characteristic phase in multiple scattering, due to the 

Faraday effect in the effective medium accumulated over a distance An ensemble of Rayleigh scatterers (for which the 
electromagnetic field changes only slightly on a scale comparable to the particle size) has the finite value 

'^^^ (27) 



(m2 + 2)2 



The Faraday effect of a composite material made of particles smaller than the wavelength, and of different shape (spherical, 
needle like, plate like) was discussed in Ref. p^ , using a more general version of the effective medium approximation (not limited 
to dilute samples). 



For Rayleigh scatterers (a; ^ 0), one obtains 



do = 



X = 



4(m2 - 1)2x3' 

27£2^ 



so that, by definition (|l7|), 



8(m2 - lYx 
Veff _ 9\/2to 



fVo (m2 + 2)2 



(28) 



Apart from a factor depending on the index of refraction, the effective Verdet constant is found to be the product of the 
volume fraction of the particles by their Verdet constant. As noticed before 0, one finds a factor of \/2 more than expected on 
the basis of the effective medium approach of Eq. ( p7|) if Veff would have defined by Sa = 2VeffB/uj. This discrepancy is due 



to the denominator in Eq. (21), 1 — (cosS)^ = 1/2 for Rayleigh particles. 

From the experimental parameters described in the experiments of Erbacher et al. (a relative index of refraction of m = 1.15), 
the estimate using Eq. (|8|) is 

« 1.32. (29) 

This value is to be compared with the experimental value of 1.55±0.15 Q]. The proximity of the two values probably explains the 
success of previous theories based on Rayleigh scatterers, although this experiment dealt with Mie particles. In this experiment, 
the maximum of the distribution of the size parameters was roughly estimated at a; ~ 23 but the width of the distribution was 
very broad. Using the parameters of Erbacher et al., our Mie theory reproduces the limit of Eq. ( p9| ) for a: = 0, but predicts a 
value of only 0.4 at a; ~ 23. The solid line in figure ^a represents the effective Verdet constant as a function of the size parameter 
X, for the same relative index of refraction used in the experiment of Erbacher et al.. One can clearly see in this figure that a 
distribution of large spheres of size parameter x of the order of 20 or higher can not explain the experimental result. However if 
the size distribution of the scatterers was rather centered around a size parameter of roughly 10, the experimental value of the 
effective Verdet constant could possibly be recovered from this Mie theory. In figure ^3, an higher value of the relative index 
was chosen. In this case, the effective Verdet constant is seen to be enhanced by resonances, which will be the subject of the 
section ^ 



6 Faraday rotation for multiple Mie scattering 

The first term in Eq. (|2l|) originates from the Faraday effect inside the scatterers and is the main contribution in the Mie regime 
of a; ^ 1. Resonances will be discussed in section 0. Using the definition of Eq. (p7|), the effective Verdet constant of Mie 
particles is written as 
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Veff _ SmQscattVSX , , 

When the size parameter obeys x S> 1 and y — mx 3> 1, the scattering can be interpreted in terms of geometrical optics. 
In geometrical optics, rays incident on the sphere are considered rather than plane waves. The ray with central impact is 
characterized in Mie theory by J = 1. In figure ||a, the effective Verdet constant, that contains contributions from all the rays, is 
plotted as a solid line, with respect to the size parameter x. The dashed line in this figure represents the separated contribution 
of the ray with central impact J = 1 in the effective Verdet constant. The two curves merge for x — 0, the Rayleigh limit, but 
deviate from each other for larger size parameter. For a given value of the size parameter, it can be noted that the ray with 
central impact already represents a significant contribution to the effective Verdet constant. In addition to this, the effective 
Verdet constant is seen to decrease for increasing size parameter x. This is an important observation, because this means that 
Mie scatterers of large size (typically x > 20) are less efficient than Rayleigh scatterers in suppressing the coherent backscattering 
cone. 



7 Mie resonances 

Let us first recall some results for resonant Rayleigh scattering, for which the resonance behavior is analogous to the resonance 
of a two-level atom in atomic physics. For resonant Rayleigh scatterers the effective Verdet constant is related to the "path 
length" Lpath of the wave inside the particle. In fact, the path length, the dwell time of the light and the total electromagnetic 
energy stored inside the scatterer have been seen to be proportional. The path length can be defined as (see formula (2-22) P17 
ofRef. (l5)), 

Lpath = lim (31) 

where Qabs is the absorption cross-section and rrii the imaginary part of the index of refraction. The physical idea behind this 
definition is that the longer the path of the light is in the particle, the more the light will suffer from absorption. For resonant 
Rayleigh particles, the relation to the effective Verdet constant is obtained from Eqs. ( p8pl| ), 

fVo 3 a ' 

At resonance, the path length can exceed the size of the scatterer, which means that the effective Verdet constant should be 
strongly enhanced by resonant scattering. Alternatively, one may relate the path length to the time spent by the wave in the 
medium, which means that the Faraday rotation is in some sense a Larmor clock, measuring this time 

The question is whether resonant enhancement of Faraday rotation occurs in resonant Mie scattering. For Mie resonances 
the increase in path length is well related to the change of the electromagnetic energy within the scatterers with respect to the 
surrounding. The total time-averaged electromagnetic energy inside the sphere is denoted by W and Wo represents this energy 
for the incident plane wave. For weak absorption, mi <C Wr, the electromagnetic energy W can be approximated by [l^ 

W ^ SrUrQabs _ SrUrLpath ^^^^ 

Wo Sxrui a 

where Eq. ( |3l| ) was applied to obtain the last equality. This relation is exact for scalar waves and is a very good approximation 
for vector waves. It is even an excellent approximation in the vicinity of resonances were the deviations between the exact 
solution and its approximation are the largest. In the domain of J ~ a;, several resonances take place in the Mie coefficients aj 
and bj. These resonances are well separated and can be numbered by an additional integer k the order of the resonance p^ . 
Near one electric Mie resonance of a specified order, the denominators of the Mie coefficients aj and cj, which are identical, are 
close to zero. From Eqs. (21 2^), one finds a/A ~ |cjp/|aj|. By Eqs. ( ^ and (p3|), this implies 



^H.,n.l~^,<.,„^^. (3.) 

In the last equality the role of \aj\ is dominated by |cjp at resonance, since resonances in the scattering cross-section Qscatt are 
much less significant than resonances in W. Indeed, in Fig. (^), the lower curve represents Qscatt , on normal scale, near one 
resonance of a water sphere. It is much below the curves of W or Ve///(/Vb) which are even plotted on a logarithmic scale. From 
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Eq. (^), the effective Verdet constant is expected to be simply proportional to the electromagnetic energy W (or equivalently 
to Lpath)- This generalizes the result for Rayleigh scatterers in Eq. (|3^). The numerical verification of Eq. ( ^4[ ) can be deduced 
from the double logarithmic plot of ^///(/Vb) against the total electromagnetic energy W in Fig. (^). 

The proportionality of the Verdet constant and the total electromagnetic energy (or equivalently Lpath) has been derived for 
the particular case where the path of the light is confined along the same line (ID problem as in a Fabry-Perot configuration 
for instance |^). In this case, the cumulative character of the Faraday rotation with respect to the path length leads to an 
experimentally observed enhancement of the Faraday rotation [p^ . Eq. ( 3^) applies to any resonant impact and shows that the 
possible occurrence of spin flips in Mie scattering, as suggested in Ref. \7], does not affect the behavior of the effective Verdet 
constant near resonances and that the Faraday rotation still accumulates along the path as in the ID case. This is consistent 
with the observation made in section |[ that the ray with central impact had an important role for interpreting the Faraday 
rotation for multiple Mie scattering. 

In conclusion, like for resonant Rayleigh particles, a strong correlation between the effective Verdet constant and the stored 
energy inside the sphere was found for resonant Mie particles. The general behavior and proportionality in the vicinity of a 
resonance is apparently universal. 

8 Shift of the intensity profile of the coherent backscattering cone in a mag- 
netic field 

A magnetic field can be expected to introduce some anisotropy in the light intensity profile of the cone. In this section, the 
form taken by this anisotropy is investigated by taking care of the selection rules imposed for the polarization in reflection of a 
semi-infinite system of Mie scatterers. This analysis is restricted to linear corrections in the magnetic field, so that the parameter 
A discussed in the former section, quadratic in the field, will no longer appear. 

This approach is based on an improved version of the scalar diffusion approximation. The ladder propagator at point 
r — {r±,z} for a source at r' = {r'j^,z'} in a semi-infinite medium is denoted by p(r,r'). The z-axis is directed along the 
normal of the sample, and rj^ and r'^ are vectors perpendicular to the z-axis. Because of translational invariance in the plane of 
the sample, the ladder propagator only depends on rj_ — r'^. The two-dimensional Fourier transform of p(r,r') with respect to 

— is denoted by p(q, z, z'). The ladder propagator p(r, r') obeys the following diffusion equation 

(-V2 + -l)p(r,r')=<5(r-r'), (35) 



with the radiative boundary condition 



Vrx,Vrl,Vz>0, p({ri, z}, K, z' = -zq}) - 0. (36) 

The trapping plane is located at a distance zq = 2€*/3 outside the sample, and La is the absorption length for the light 
intensity. We assume that the first and the last scattering events take place one transport mean free path i* away from the 
boundary, in the directions specified by the incoming wavevector p and the outgoing wavevector p'. This allows to calculate 
the contribution of the ladder diagrams Lj-^;(p, p', q, B), the so-called incoherent contribution to the coherent backscattering, 
directly from p(q = 0, z — £*pz,z' + (*p'^). Similarly, the coherent contribution is obtained from the most-crossed diagrams 
Cj^;j,(p, p', q, B), which are derived from L~j.;(p, p', q, B) using Eq. (p^. These derivations are detailed out in Appendix A. 
The coherent contribution in an helicity basis and at reflection (z = 0, z' = 0) reads, 

C+,> =p{Y>^+v'^,-tp,,t'p'^) X 5„„, X (l + 2feicoMet(pl,pi,B)-2(762^^^B-(pi+pl)). (37) 

where the Kronecker symbol 5aa' guarantees conservation of helicity cr(p) = cr'(p'), and 6i and 62 are real-valued coefficients 
to be determined. Most experiments on coherent backscattering were done in the helicity-conserving channel, which has the 
advantage of having a maximal enhancement factor (since the contribution from single scattering vanishes in this case) and of 
having an isotropic line shape. Eq. (^7|) states that the magnetic field modifies the cone exactly in this channel, in agreement 
with previous work Even when a magnetic field is present, no coherent backscattering is found in the opposite-helicity 
channel, at least according to the present diffusion approximation. 

Only the components of the magnetic field along the slab contribute in the r.h.s of Eq. (|37|). When the field is perpendicular 
to the slab, the decrease of the enhancement factor, described by A, is the sole impact of the magnetic field and is not included 
in the present approximation. When the magnetic field is in the plane of the slab, two corrections show up. The first one, 
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proportional to 61, is magneto-transverse, since it produces a shift of the intensity profile of the cone in the plane of the slab, 
normal to the magnetic field. This correction is independent of the state of helicity of the light. The second correction in Eq. 
( |37| ) proportional to &2 does depend on the helicity a. It produces a shift of the intensity profile of the cone in the direction of 
the magnetic field in the plane of the slab, quite similar to the correction induced by the magnetic field in the group velocity 

. . . n 

The coefficient 61 can be calculated independently from 62, in a way exactly analogous to the calculation of ai in Ref. 
responsible for the Photonic Hall Effect (PHE). The result is 



1 /!i^cosgsingE..>5m(r,V(^)ri^,(B^g,. ^ 

{^-{^o,e)f 2/!^dcos0E..' \TL.m' ' ^^^^ 



where the magneto T- matrix Tl^,{R = g)(0) was introduced in Eq. (pi)[). Note that the imaginary part in Eq. ( p8| ) discriminates 
hi from the parameter ai responsible for the PHE, in which the real part figures. For Rayleigh scattering, one can readily prove 
that the parameters 61 and 62 both vanish. The calculation of 62 for Mie scattering is very complicated and is beyond the scope 
of this paper. 

The expected modification of the lines of equal enhancement factor due to the magnetic field is now investigated. Equation 
(I37I ) translates into a CBS line shape 

E{ii, (f) = 1 + C(;u)(l + 2bifismip — 2b2cr fi cos tp) . (39) 

The enhancement factor, in the helicity-conserving channel, is denoted by E for a state of helicity a. The dimensionless parameter 
H = (jj£0 and C (/i) the well-documented line shape of the cone without applied magnetic field were introduced Q . The azimuthal 
angle between the projection of the outgoing wave vector into the plane of the slab and the magnetic field direction, which has 
been chosen along the x-axis, was denoted by (p. For simplicity, only the magneto-transverse correction proportional to bi will 
be considered here, so that 62 = 0. This case corresponds to unpolarized incident light, for which the term proportional to 62 
in Eq. ( |39| ) vanishes. The pattern of the lines of equal enhancement factor associated with the &2-correction alone is the same 
as the one of the transverse correction after a rotation of angle tt/2 about the x-axis. For a typical experiment |^ with CeFs 
particles of approximate radius of 2/im at room temperature, with a wavelength A = 457nm, Eq. ( |3^ ) leads to the estimation 
61 ~ 1.8 • 10~^. The experimentally measured mean free path for a volume fraction of / = 0.1 is i* ~ 90/zm. Eq. ( |39| ) is valid for 
<C 1, which means that our approach is limited to the angular domain \9\ < 0.8rad. The equation for the lines of constant 
enhancement factor Eq in the absence of magnetic field is i?o = 1 + C{ii^), independent of the azimuthal angle Lp. The first order 
correction in the magnetic field is separated by writing /i = /J^ + fi^. Equation (^9|) gives 

M' = -2foig^/sin^. (40) 

In Fig. (||), the polar diagram of the lines of constant enhancement factor is shown for 61 = 1.8 • 10^^. As apparent from Eqs. 
( |39| , ^ ), the distortion of the lines should increase away from backscattering (at exact backscattering there is no modification at 
all to first order in the magnetic field). As a consequence, with the value of bi given above only a modification of the line shape 
in the wings might be observed. Although the condition bifi <^ 1 limits the domain of validity of the approximation there should 
nevertheless be a sufficiently broad angular range, in the wings of the cone, where the magnetic corrections could be visible. In 
figure 1^ the condition of validity of the approximation has been satisfied. 

Another condition of validity lies in the use of the diffusion approximation. This approximation predicts a l//x^ behavior 
for the line shape of the cone in the wings, which is actually a wrong result. In the wings, the contribution of lowest orders of 
scattering is dominant and not properly taken into account in the diffusion approximation. Using an exact theory, Gorodnichev 
proved that the outcome is a 1/^ dependence His result has been derived only for point-like scatterers but it should also 
be valid for Mie scatterers when the mean free path is much larger than the wavelength. In any case, the magnetic correction in 
Eq. (^) depends only on the logarithmic derivative of C^jjP), which should change only by a factor of two if the actual law is 
l//i or l//i^, the general pattern of the lines of equal enhancement factor being not modified. Therefore, the shift of the center 
of mass of the light intensity profile which was calculated, should be fairly robust with respect to the exact form of C{f/*). 

9 Conclusion 

This paper describes two modifications on the coherent backscattering cone produced by a magnetic field. The first one, the 
decrease of the enhancement factor, depends on the parameter A quadratic in the magnetic field, and was observed experimentally. 
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The second modification is related to the anisotropy of the light intensity and appears already in linear order of the magnetic 
field. Preliminary experiments seem to have reported the possibility of a shift in the intensity profile [2^ . Our analysis applies 
to spherical scatterers of any size that are Faraday active. The decrease of the backscattering cone gets less pronounced as 
the size of the scatterers increases, whereas the shift in the intensity profile is only possible with finite size scatterers. As was 
surmised in Ref. Q, the effective Verdet constant defined from the decrease of the cone is enhanced near Mie resonances. The 
effective Verdet constant is found to be intimately related to the stored electromagnetic energy i. e the dwell time of the light in 
the particle. 

We acknowledge R. Lenke for making available his recent experimental work on the effect of a magnetic field on coherent 
backscattering light. We thank JJ. GrefFet for stimulating discussions. 



A Derivation of the magnetic corrections to the tensor C^^, 



In this appendix, the notations are explained and a demonstration of Eq. (|3^) is given. The transverse part of the free Green 
tensor is denoted ^ 

Go = 7 — 7 — NO . (41) 

with (Ap)y = 5ij —piPj the projector upon the space orthogonal to p. Similarly, the hermitian projector on the space orthogonal 
to p for a given state of helicity a is -Pj'^(p) — 5 ((Ap)^ — i<J£ijkPk) ■ 

Generalizing Eq. (|lj) for finite value of q (to first order in q and in the magnetic field B) gives 

d(p, q, B) = I + + [L(p, q) - r- (p, q, B)] , (42) 

IUjJ 

where L(p,q) = 2(p • q)I — pq — qp and r~(p, q, B) are tensors of rank two, linear in q, that determine the anisotropy in 
diffuse scattering. Without a magnetic field it is well known that r"(p, q) = 2(p • q)/(l — {cosff)). When a magnetic field is 
present, the first order correction in the magnetic field is separated as r~(p,q, B) = r'^(p,q) + _B(5r^(p, q, B). Because of 
the symmetry relation of Eq. (^), one has r~(p,q,B) = r~(p,q,— B)*. This implies that (5r^(p,q,B) = — (p, q, B)*, 
so that (5E~(p, q, B) must be anti-hermitian. Mirror symmetry imposes in addition that Tpp'(B) = T_p_p'(B) and thus 
E^(— p, — q, B) = E^(p, q, B). The general form of the tensor 5E+(p, q, B) allowed by symmetry has already been discussed in 
Ref. ||l^ and will be only slightly different for 5V^: 

(5Fj^ (p, q, B) =z6idet(p,q, B)(5ij -ieijfcPfc 162(6 • q) + ^63(6 • p)(p • q) + 164 [pfc^'feigj + giPfe^fcj] , (43) 

where the bi are real-valued coefficients to be determined. 

The contribution of the first and the last scattering events is obtained by multiplying the tensor for the ladder diagrams of 
Eq. jl^ ) by free Green tensors Go. This gives the incoherent contribution to the coherent backscattering, 



GoGS-L+(p,p',q)-GoGS ~ Go-d(p,q)-GS G*.d(-p',-q).Go, 

~ (1 - p • q)p(q = 0, z, z') (1 + p' • q) Ap Ap, , 

~ p(q = o,z-rp„z' + rpi)ApAp,, 

in terms of the Fourier transform of the ladder propagator p(q, z,z'), which is obtained from the solution of the diffusion 
equation (|35|). The coherent contribution depends on the most-crossed diagrams 0^^^. The tensor C^u^ is obtained from i.^j,; 
after reversing the indices k and Z, according to Eq. (|l^), and after adding the contribution of the first and the last scattering 
event. In this calculation of the tensor C^^j^, the reciprocity transformation ( p^ for the normal components z and z' of the 
tensor has been neglected. The result is finally evaluated at reflection where z = z' = 0. The components of the ingoing and 
outgoing wave vectors p, p' perpendicular to the z axis are denoted by pj^, p'^. Finally the coherent contribution can be written 
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AiiAjk + i 



jk 



-Pa 



P±r 



B 



Mr, 



-,P± + p'l,B Akj 



with the definition 



M^(p,q,B) = -C$, 



(p,q,B), 



(44) 



(45) 



and C = 2(irf2 ^ v) is the Mie generahzation of the parameter F in Eq. (73) of Ref. |Q, which was shown to produce a rotation 
of the polarization vector in the hnear polarization channels of coherent backscattering. To first order in the magnetic field, the 
front factor p(p_L + p^, —^*Pz, i*p'z) of Eq. fi4) is evaluated by replacing the parameter La, which was an absorption length for 
the evaluation of the light intensity in Eq. (|3q), by a factor depending on the backscattering angle according to: 



1 



(PJ 



In second order in the magnetic field, A would be present here as well as stated in Eq. (|T^). 

The values of the coefficients bi can only be found by solving a system of four coupled equations that one obtains when 
inserting Eqs. ( ]4^ , ^ ) into Eq. (^9|), and which is not reported explicitly here. The contribution of 63 always vanishes in Eq. 
(Q), since it is proportional to the scalar product p • q in Eq. (^3[), which is transformed into (p' — p) • (p + p') = in the 
operation involved in Eq. (|T^). For the same reason there is no contribution of r'^(p, q) in Eq. (Q). Selection rules for the 
polarization are obtained in the helicity basis by considering the product 



prkip)cLiPuip')- 



(46) 



In this calculation, the terms proportional to C in Eq. (|4g) disappear, as well as the contribution from 64 which is longitudinal 
as can be seen from Eq. (|4^). Among the four terms of (5r(p, q, B)^, only the terms proportional to bi and 62 survive, and Eq. 
Wn) is obtained for Cj\, . 
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i,p + q/2 X X X X X X j,p'+q/2 



k,p- q/2 X X X X X x l,p'- q/2 

+/- B 

+ B 

i,p + q/2 X— X- — X X — -X— X j,p'+q/2 



k,p- q/2 X x^ X X X x l,p'- q/2 

+/- B 



+/- 



= Lijki (p,p'lq,B) 



+/- 

= Cijki (p,p'lq,B) 



Figure 1: Ladder diagrams Lf^j^i and the most-crossed diagrams Cf^if. in a magnetic field. Bold lines denote the Dyson-Green 
tensor, the crosses denote T-matrices and dotted lines connect identical particles. 
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Figure 2: Asymmetry parameters (cos0) and (cos( 
m = 1.33 (a) and m = 2.73 (b). 
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Figure 3: Plot of the effective Verdet constant as a function of the size parameter x. (a) The solid line represents the 
solution of Eqs. (pc|-p3|), containing contributions from all the rays incident on the sphere, whereas the dashed line corresponds 
to the contribution of only the first partial wave J — 1 (the ray with central impact in geometrical optics). The scattering 
medium has an index of refraction of 1.7, the value in the experiment of Erbacher et ai. (b) The same plot for an index of 
refraction of 2.73 for which resonances are clearly visible. The host medium is in both cases glycerol (of index of refraction 1.47). 
A general decrease of the effective Verdet constant for increasing size parameter x can be observed in both plots. 



13 



(a) 



(b) 




47.060 47.065 47.070 47.075 0.3 0.8 1.4 1.9 

Size parameter x log(WAV„) 

Figure 4: Near a particular resonance curve of a water sphere of index m = 1.334 — 1.5 • 10~^i the following curves have been 
plotted: (a) Respectively from the upper part of the figure to the bottom: a plot of log(ye///(/Vb)) (solid), logW^ (dashed), 
and the scattering cross-section Qscatt (dotted) as a function of the size parameter x. (b) Close to this particular resonance, 
log(Ve///(/Vb)) is shown against the total electromagnetic energy W (dots), to be compared with a line of slope one (solid), the 
prediction of Eq. (^J). 
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Figure 5: Polar plot of the lines of constant enhancement factor of the coherent backscattering cone in a magnetic field, in the 
helicity- conserving channel, for bi = 1.8 • 10~^. The magnetic field is along the horizontal axis on this graph. The graduations 
of the axis represent the dimensionless parameter /i = uj(.6. For instance, a graduation of 10 corresponds to an angle of 145mrad 
(for the mean free path I* ~ 90/im, and the wavelength A = 457nm mentioned in the text). Without magnetic field or for a 
magnetic field perpendicular to the slab, these lines of constant enhancement factor would have been circles. 
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